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A method of and apparatus for processing seismic data 

The present invention relates to a method of processing seismic data. It particularly 
relates to a method of processing seismic data that includes two corresponding modes of 
seismic energy propagation, for example where one mode arises from partial mode 
conversion of another mode at a boundary of a layer of the earth. The method of the 
invention provides an estimate of static shifts caused by a layer at or near the earth's 
surface that has different seismic properties from the underlying layers. The invention 
also relates to a method of seismic surveying that includes processing acquired seismic 
data in the above way. The invention fiirther relates to an ^paratus for processing 
seismic data. 

Figure 1 is a schematic view of a seismic surveying arrangement. In this figure the 
surveying arrangement is a marine surveying arrangement in which seismic energy is 
emitted by a seismic source 1 that is suspended beneath the sea-surface from a towing 
vessel 2, When the seismic source 1 is actuated seisrnic energy^is emitted downwards, 
• passes into the earth's interior and is reflected by a geological feature 3 that acts as a 
partial reflector of seismic energy. The reflected seismic energy passes upwards 
through the earth's interior, into the sea and is detected by an array of seismic receivers 
4 disposed on the sea-bed. The seismic energy path shown in Figure 1 is generaUy 
referred to as the '^prinaary path". Information about the earth's interior is obtained, for 
instance, by determining the travel time of seismic energy from the seismic source 1 to 
the receiver 4 along the primary path. 

In practice refraction wiU occur at the seabed and at boundaries between layers of 
different geological composition within the earth, but this has been omitted from Figure 
1 for clarity. Only one reflector of seismic energy is shown in Figure 1 but, in practice, 
a number of geological features within the earth's interior will act as partial reflectors 
for the seismic energy. Moreover, only one source 1 is shown in Figure 1 but m 
practice a marine seismic surveying arrangement may have an array of sources. 
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One problem encountered in processing seismic data is that the geological structure of 
the earfli is not uniform. In particular the properties of the region 5 of the earth near the 
surface (hereinafter referred to as the "near-surface") are j&equently significantly 
different from the properties of tihe underlying geological structure 6 (hereinafter 
referred to as the "basement"). This can occur if, for example, the near-surface is less 
consolidated than the basement In particular, the velocity of seismic energy may be 
significantly lower in the near-surface 5 than in the basement 6, and for this reason the 
near-surface is sometimes known as a "low- velocity layer" (LVL). This difference in 
velocity will produce a shift in the travel time of seismic energy from the soxirce 1 to a 
receiver 4 compared to the travel time that would be recorded if the near-surface and the 
basement had identical seismic properties, and these shifts in travel time are generally 
known as "static shifts", or just "statics". 

The near surface may also aJBfect tihe seismic energy incident on the receiver in other 
ways. For example, the near surface may have a focusing or defocusing effect on 
seismic energy propagating through it, and this can affect the amplitude of seismic 
energy received at a receiver. The near-surface may also attenuate the amplitude of 
seismic energy propagating through it. Furthermore, anelastic attenuation and 
focusing/defocusing effects in the near surface may also distort the waveform of the 
seismic energy wavelet. 

Although Figure 1 illustrates the near surface 5 in connection with a marine seismic 
survey, a similar near-surface region occurs on land so that static shifts also arise in a 
land-based seismic survey. 

The near surface 5 is shown as extending to the earth's surface (the sea-bed) in Figure 1, 
but it need not extend to the surface and there could be a fiirther layer (not shown) 
overljdng the near surface 5. Furthermore, although Figure 1 shows a well-defined 
boundary between the near-surface 5 and the basement 6 in practice there may be a 
gradual transition between the basement 6 and the near-surface 5 rather than an abnq)t 
boundary. 
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The static shift generated by the near-surface 5 depends on the thickness of the near- 
surface, and on the velocity of propagation of seismic energy through the near-surfece. 
Lateral variations usually occur in both the thickness of the near-surface 5 and the 
propagation velocity through the near-surface, so that the static shift observed at a 
seismic receiver at one location is likely to be different ftom the static shift observed at 
a receiver at another location. To a first approximation, the entire data set recorded at 
one receiver will be advanced or delayed by a static time shift relative to data recorded 
at another receiver. 

It is highly desirable to take account of the static shift when processing seismic data. 
Unless these static shifts are removed from the seismic data, ambiguity will exist as to 
whether variations in airival times of seismic events from deeper layers are due to 
variations in the depth or lateral locations of those deeper layers, or simply arise owing 
to propagation effects in the near-surface 5. 

It will be seen from Figure 1 that the seismic energy passes through the near-surface 5 
twice, once as it travels from the seismic source 1 to the reflector 3, and again as it 
travels from the reflector 3 to the receiver 4. Both traverses of the near-surface 5 wiU 
cause static shifts; the shift caused as the downwardly propagating seismic energy 
passes through the near-surfece 5 is known as the source-side static shift, and the shift 
caused as the upwardly propagating energy passes through the near-surface 5 is known 
as the receiver-side static shift. 

The method of the present invention makes use of an s-mode event and a corresponding 
p-mode event to obtain information about flie low velocity layer. It assumes that liie 
differences between tiie s-mode event and the corresponding s-mode event, for example 
differences in the travel time, ampUhide and/or waveform of the events, arise wholly or 
primarily in the near-surface 5. 

One embodiment of the invention makes use of an s-mode event and a corresponding p- 
mode event where one of the modes has been generated by partial mode conversion of 
the other mode. For example, if tbe contrast in physical properties between tiie . 
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basement 6 and the near-surface 5 is sufiSciently great, signijBcant conversion between a 
p-mode of energy propagation and an s-mode of energy propagation, or vice-versa, 
takes place as seismic energy propagates through the boundary between the basement 6 
and the near-surface 5. This mode conversion is described more fully below with 
reference to Figure 2. As a result of this mode conversion, for every p-mode event 
recorded by a seismic receiver above the boundary between the near-surface and the 
basement there will be a corresponding s-mode event. The p-mode event and the 
corresponding s-mode event will occur at a receiver 4 at different times after actuation 
of the seismic source, because the p-mode and s-mode propagation velocities ia the 
near-surface are not equal to one-another. However, since the mode conversion occmrs 
at the lower boundary of the near-surface, the time delay between the p-mode event and 
the corresponding s-mode event must arise as the result of different velocities for the 
two modes in the near-surface. The time delay between the p-mode event and the 
corresponding s-mode event will not be significantly influenced by the properties of the 
basement. The travel time difference between the p- and s-modes through the near- 
surface can thus be readily detenrdbtied by identifying a p-mode event and the 
corresponding s-mode event in the seismic data and determining the time delay between 
the p-mode event and the corresponding s-mode event. 

For any parameter indicative of an aspect of the seismic data such as, for example, the 
pressure or a component of the particle motion (such as a component of the particle 
displacement or particle velocity), events of one of the p- and s-modes will generally 
appear more strongly than will events of the other mode. Thus, in general, a p-mode 
will occur predominantly in one parameter, such as the vertical component of the 
particle motion and the corresponding mode-converted s-mode will occur 
predominantly in another parameter, sucTi as a horizontal component of the particle 
motion. 



Mode conversion may also occur when a downwardly propagating wave undergoes 
reflection at the interface between the near-surface S and the basement 6. The iavention 
may be applied to this case since, where mode conversion occurs on reflection, the 
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reflected signal will contain a p-component and a corresponding s-component, and the 
two components wiU have different travel times through the near-surface 5. 

Mode conversion may also occur when a downwardly propagating wave undergoes 
critical refraction at the interface between the near-surface 5 and the basement 6, to 
generate a seismic wave propagating along the interface between the near-surface 5 and 
the basement 6. The wave propagating along the interface will excite both p-mode and 
• s-mode waves in the near-surface, and the invention can be apphed to these p- and s- 
modes. 

Co-pending UK patent appHcationNos. 0101936.3 and 0110573.3 disclose a method of 
determining the difference betweai the static_d^^ and the static shift for 

an s-wave generated by mode-conversion at the boundary between the near surface and 
the basement. Pairs of a p-event and a corresponding s-event are identified in the 
seismic data, and the travel time difference between the p-event and the corresponding 
. s-event is found,. This travel time difference.is the difference between the s-wave static 
shift and the p-wave static shift for a particular receive:. 

Many seismic surveys now use multi-component receivers that record two or more 
component of the seismic energy incident on the receiver. For example a 3-component 
(3-C) seismic receiver can record the x-, y- and z-components of the particle motion at 
the receiver, and in a marine survey a 4-component (4-C) seismic receiver can record 
the pressure of the water column (which is a scalar quantity) in addition to the x-, y- and 
z-components of the particle motion. A further problem arises in processing multi- 
. component seismic data, which is that a receiver may record some components of the 
seismic data more accurately than other components. This problem is known as vector 
infidehty. One cause of vector infideUty is inadequate couphng between the receiver 
and the earth's surface so that, m a receiver array, the effects of vector infidelity may 
vary from one receiver of the array to another. Such vector infidehty will degrade the 
rjssults obtained by processing the seismic data. 
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A first aspect of the present invention provides a method of processing seismic data 
including corresponding first and second modes of seismic energy, the method 
comprising the steps of: 

a) . processing a trace indicative of a first parameter of the seismic data acquired at a 
first receiver including at least events corresponding to the first mode with a trace 
indicative of a second parameter of the seismic data acquired at the first receiver 
including at least events corresponding to the second mode; 

b) identifying an event in the processed data corresponding to partial mode 
conversion; and 

c) obtaining information firom the amplitude and/or the waveform in the processed 
data of the event corresponding to partial mode conversion. 

The present invention may be used to obtain information about properties of the near 
surface firom the amplitude and/or the waveform of the event corresponding to partial 
mode conversion. Additionally or alternatively, it may be used to obtain information 
about variations, in coupling strength &om one receiver to another. 

Iq a preferred embodiment the invention is applied to data in which the second mode 
was generated by partial mode conversion of the first mode at a boimdary face of a layer 
of the seabed. However, the invention is not limited to this, and may generally be 
applied to any pair of events, one in one mode and one in another mode, for which 
differences between the amplitude and waveform of the two events will arise primarily 
in the low velocity layer. 

The parameters of the seismic data that are used in the invention are preferably chosen 
so that a p-mode will occur predominantly in one parameter and the corresponding 
mode-converted s-mode will occur predominantly in the other parameter, 'riie 
parameters may be parameters that are acquired directly at the receiver, such as the 
vertical component of the particle motion and a horizontal component of the particle 
motion. Alternatively, the parameters may be obtained by processing the seismic data 
acquired at the receiver, for example by decomposing the acquired seismic data into its 
up-going and down-going constituents or into its p- and s-constituents. 
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The method of the inventioii may be used to process pre-existing seismic data. 
Alternatively, it may be incoiporated in a method of seismic surveying for processing 
the data as it is acquired or subsequently. 

A second aspect of the present invention provides a method of seismic surveying 
comprising the steps of: directing seismic energy propagating in a first mode towards a 
boundary face of a layer of the seabed such that partial mode conversion of the seismic 
energy into a second mode occurs at flie boundary face; acquiring seismic data 
including the first and second modes of seismic energy at one or more receivers; and 
processing the seismic data according to a method as defined above. 

A third aspect of the invention provides an apparatus for processing seismic data 
including first and second modes of seismic euergy, the apparatus comprising: 

a) means for processing a trace indicative of a first parameter of the seismic data 
acquired aLa.firetjeceiyerjncluding at lejst events correspondm^ 

a trace indicative of a second parameter of the seismic data acquired at the first receiver 
including at least events corresponding to the second mode; 

b) means for identifying an event in the processed data corresponding to partial 
mode conversion; and 

c) . means for obtaining information from the ampUtude and/or the waveform of the 
event in the processed data corresponding to partial mode conversion. 

The apparatus may comprise a programmable data processor. 

A fourth aspect of the present invention provides a storage medium containing a 
program for a data processor of an apparatus as defined above. 

Preferred embodiments of the present invention will now be described in detail by way 
of illustrative example with reference to the accompanying figures in which: 



Figure 1 is a schematic view of a marine seisanic survey; 
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Figure 2 is a schematic illustration of mode conversion of an up-going p-wave at the 
boundary between the basement and a low-velocity layer; 

Figure 3 shows two schematic examples of seismic data traces where the mode 
conversion of Figure 2 occms; 

Figure 4 shows the result of cross-correlating the two seismic data traces of Figure 3 
with one another; 

Figure 5 illustrates results of the cross-correlation method of UK patent application Nos. 
0101936.3 and 0110573.3; 

Figure 6(a) is an enlarged view of selected traces from the results of Figure 5; 

Figure 6(b) shows, the results, of normaUsing the.traces_.of Rgur 6(a); 

Figure 7 shows the variation in amphtude of a selected event from the traces of Figure 
6(b); 

Figure 8 is a block flow diagram of a method according to one embodiment of the 
present invention; and 

Figure 9 is a schematic block diagram of an apparatus according to one embodiment of - 
the present invention. 

The principle of the invention will now be described with reference to an embodiment 
in which the pair of corresponding events are generated by partial mode conversion. 
Figure 2 illustrates the occurrence of mode conversion as an upwardly propagating p- 
wave passes through the interface between tiie basement 6 and the near-surface 5. The 
up-going p-wave 7 is partially transmitted as a p-wave 7', but also undergoes partial 
mode conversion to an s-wave 7". It is assxmied that the upward-going p-wave was 
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generated by a reflection from a deep interface between two layers of the earth's 
interior, or by refraction through deeper layers. Refraction will also occur at the 
interface between the basement 6 and the near-surface 5, and the angle of refraction for 
the transmitted p-wave T will be different from the angle of refraction for the converted 
s-wave 7". 

The transmitted p-wave 7' and the converted s-wave 7" will both produce a signal at the 
receiver array. The transmitted p-wave T and the converted s-wave 7" are. as indicated 
in Figure 2, travellmg in different directions and hence will arrive at different receiver 
positions,notatthesamereceiver. However, a similar converted s-wave from a 
neighbouring conversion point will arrive at the same receiver as the transmitted p-wave 
T shown in Figure 2. It will be assumed that, for a transmitted p-wave and a converted 
s-wave received at a particular receiver, the difference between Ihe conversion point (at 
the base of the near-surface) of the converted s-wave and the point in the base of the 
near-surface through which the p-wave was transmitted does not significantly affect 
estimates.ofthe.staticshifts.obtained>yltieme*odpfJM J^^^^ 
acquired at the receiver will thus contain an event arising from the arrival of the 
transmitted p-wave 7' and an event arising from the arrival of the corresponding mode- 
converted s-wave 7". 

The transmitted p-wave 7 will be recorded predominantly on the vertical component of 
particle velocity recorded at the receiver, while the converted s-wave 7" will be 
recorded predominantly on the radial component of the particle velocity recorded at the 
receiver 4. (The radial direction is the source-receiver direction projected onto the sea- 
bed, and this direction wiU be defined to be the x-direction.) This is because in a 
practical seismic surveying arrangement waves that have propagated up from deep in 
the earth's interior tend to make an angle of 30° or less to the vertical for a typical 
structure of the earth's interior and a typical depth of the target reflector. The shear 
wave is recorded predominantly on the radial con5)onent because the particle motion for 
a shear wave is perpendicular to the direction of propagation, compared to the particle 
motion for a p-wave which is along the direction of propagation. 



M&C Folio No P51798GB " 10 * 

Figure 3 is a schenxatic illustration of traces indicative of two parameters of seismic data 
recorded at a receiver 4 where partial mode-conversion as shown in Figure 2 occurs. 
The two parameters shown in Figure 3 are the radial component (x-component) of the 

T^aT-Hr»1<a T/olrvo-i+rr 0+ f1l<» fooo-t-vrpT n-rtA ■fl-ic* -xrairfi r%*^1 />rxrv«-»%^<*^ A'v^-f /fy /^^<•-v-^-*-^r>»-*%#-s•*'#•^ f\-P 

the particle velocity measured at the receiver. Since the s-wave 7" was created as a 
conversion from the up-going p-wave 7, if the interface between the basement and the 
near-surface is approximately planar locally, then the recorded s-wave will contain a 
similar wavelet to the p-arrival signal. The s-wave arrival will not occur at the same 
time as the p-wave arrival but will arrive after a time delay dt (compared to the p-wave 
arrival) caused by the difference between the velocity of p-waves and the velocity of s- 
waves in the near-surface 5. This time delay dt is exactly the difference between the p- 
waye static shift and the s-wave static shift. 

In the method of co-pending UK patent application Nos. 0101936.3 and 01 10573.3 the 
first parameter is cross-correlated with or deconvolved from the second parameter. The 
resultant cross-correlograni obtained by cross-correlating the two traces of Figure 3 is 
shown in Figinre 4. There is a peak in the cross-correlogram at zero time shift:, which 
arises because the p-wave 7' arrival has a non-zero amplitude in the x-component of the 
measured particle velocity (although it occurs predominantly in the z-component of flie 
measured particle velocity). Any s-wave energy fliat appears in the vertical component 
of the measured particle velocity will also cause a peak in the cross-correlogram at zero 
time shift. The cross-correlogram also has a peak at time dt and this peak occurs from 
the p-wave arrival peak in the vertical component of the measured particle velocity and 
the corresponding, mode-converted s-wave arrival peak in the x-component of the 
measured particle velocity. Thus the time shift dt, and hence the difference between the 
p-wave static and the s-wave static, may be determined from the cross-correlogram. 

Figure 5 is a finther illustration of the results of the method of co-pending UK 
Patent Application Nos. 0101936.3 and 01 10573.3. These results were obtained 
using a receiver line containing approximately 130 receivers. Figure 5 shows, 
for each receiver, the cross-correlogram obtained from the vertical component of 
particle motion recorded at that receiver and the x-component (radial 
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component) of particle velocity recorded at thatreceiver. The difference 
between the p-wave static shift and the s-wave static shift for each receiver, as 
determmed by the method of UK Patent Apphcation Nos. 0101936.3 and 
01 10573.3, is indicated in Figure 5 as "A". The difference between the p-wave 
static shift and the s-wave static shift obtained usmg a conventional technique is 
also illustrated in Figure 5„as "B". 

. The method ofXJK Patent Apphcation Nos. 0101936.3 and 0110573.3 provides 
an unproved method of estimating the difference between the p-wave static shift 
and the s-wave static shift. However, this method makes use only of time 
inforaiation, and can detemiine only the time delay between a p-mode arrival 
and a corresponding s-mode arrival. It is not possible to determine from the 
mefliod of UK Patent Apphcation Nos. 0101936.3 and 01 10573.3 whether the 
time delay is due to the geological effect of the near-surface, or whether it is due 
to other factors such as the couphng between receivers and the earth, or other 
■sources of vector infideUty. 

The present invention sets out to derive further mformation from tiie time delay 
between a p-mode event and a corresponding s-mode event, by considering 
variations in the ampUtude and/or waveform of events in the correlation data. 

One embodunent of the present invention will now be described witii reference 
to Figures 6(a) to 7. 

Figure 6(a) is an enlarged view of selected traces fix>m Figure 5. These traces 
that show the correlation of the p-mode and the s-mode will be referred to as 
"receiver functions" for convenience. Figure 6(a) shows tiie receiver functions 
from Figure 5 for every fifth receiver, for positive values of the time lag between 
the p-arrival and the s-airival (a positive time lag indicates that the p-mode event 
arrived first at the receiver and was followed by the s-mode arrival, as shown hi 
Figure 3). 
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Figure 6(b) shows the effect of normaHsing the receiver functions of Figure 6(a), 
with respect to the amplitude of the peak in each respective receiver function at 
zero time delay. That is, each receiver function in Figure 6(b) is obtained from 
the corresponding receiver function in Figure 6(a) by normalising the receiver 
function so that the peak at zero time delay has the same amplitude for each 
receiver function in Figure 6(b). As noted above, the peak at zero time delay in 
the receiver function of Figure 4 or 5 arises through leakage of the p-wave 
arrival into the x-component of the measured particle motion, and through 
leakage of the s-wave arrival into the vertical component of the measured 
particle motion. That is, the amplitude of the peak at zero time delay in the 
receiver function for a particular receiver is indicative of the projection of the p- 
wave airiyal onto the x-component of the particle motion recorded at that 
receiver and the projection of the s-wave arrival onto the y-component of the 
particle motion. 

It may be seen by inspecting Figures 6(a) and 6(b) that the step of normalising 
the receiver functions has not had a significant effect on the amplitudes of the 
events in the receiver functions. This indicates that the amplitude of the peak at 
zero time delay is approximately the same for each receiver. This in turn 
indicates that the projection of the p-wave arrival onto the x-component of the 
acquired particle motion is approximately constant for each receiver. It also 
indicates that there is no significant variation ia the strength of the coupling for 
p-waves firom one receiver to another. 

It will be seen that the receiver functions in Figure 6(b) contain a peak at a time 
lag of around 0.3s, although the exact time lag at which this peak occurs varies 
slightly between receivers. "This pea^ the tirae delay between a p- 

wave arrival and the corresponding mode converted s-wave airival. The 
variations in the waveform and/or amphtude of this peak firom one receiver to 
another can provide information about attenuation and focusing effects of s- 
waves in the near-surface. 
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The peak in the receiver functions that arises owing to the time delay between a 
p-wave and the corresponding mode-converted s:wave will be referred to as the 
"statics pea3c" for convenience. The ampUtude and phase of the peak may also 
be influenced by variations in transmissivity and the extent to which mode 
conversion occurs at the interface between the near-surface and the basement. 
The transmissivity and conversion efficiency may vary over the interface as the 
result of, for example, changes in material properties or topography along the 
interface. 

Figure 7 illustrates the ampUtude of the "statics peak" (the peak at a time lag of 
around 0.3 seconds) for each receiver function shown in Figure 6(b), plotted 
agaiDSt.th6 recmer number. Figrae 7 show the results for every receiver 
function shown in Figure 5, and so includes data formore receivers than shown 
in Figure 6(b). 

It-will be seen-that Figure 7 exhibits considerable variation of the normahsed 
ampHtude of the peak in the receiver function from one receiver to another. 
Short-wave length variations in the ampUtude in Figure 7 (that is variations that 
occur over a small number of receivers) are Ukely to be caused by variations in 
the receiver-to-earth coupUng strength from one receiver to another. Long-wave 
length variations in the ampUtude shown m Figure 7 are more Ukely to be 
caused by geological effects arising from the near-surface structure, such as 
differences in the attenuation or focusing/defocusing caused by the near-surface 
from one receiver position to another. Tlius, analysis of the ampUtude variations 
shown in Figure 7 - for example, using a Fourier transform - can provide 
information about both geological effects and coupUng effects. 

The present invention also makes it possible to compensate for variations in the 
coupUng strength between the receiver and the earth from one receiver position 
to another. Correcting for variations in the receiver coupUiig strength is 
important, as variations in the coupling strength can distort the ampUtude and 
phase of an s-wave arrival relative to the corresponding p-wave arrival. Since 
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the s-wave arrival is recorded predominantly in the x-component of the particle 
motion required at a receiver, variations in the coupling efficiency can distort the 
recorded x-component of the particle motion, with the amount of distortion 
varying from one receiver to another. This is particularly a problem in 
^plications such as, for example, AVO ("AmpUtude versus Ofl&et) that require 
x-component particle motion recordings that are approximately consistent across 
the receiver spread and that can be directly compared to one anotiter. 

The amplitude of the normaUsed "statics peak" in the receiver functions of 
Figure 6(b) or Figure 7 is indicative of the strength of the mode-converted s- 
wave arrival which, as noted above, is caused by the near-surface. The 
amplitude of the "statics peak" is also indicative of the coupling strength of the 
receiver. The present invention makes it possible to compensate. for variations 
in the effect of the near-surface and in the coupling strength from one receiver to 
another, by normalising the amplitude of the original trace of the x-component 
of tlie particle motion acquired at a receiver relative to the amplitude of the 
"statics peak" in the normalised receiver function of Figure 6(b) for that same 
receiver. That is, for each receiver, the ampUtude of thq "statics peak" in the 
normalised receiver function is determined, and this amplitude is then used to 
normalise the original trace of the x-component of the particle motion at that 
receiver. This procedure is repeated for each receiver — the trace of the x- 
component of the particle velocity recorded at that receiver is normalised by the 
amphtude of the "statics peak" in the normalised receiver function for that 
receiver. The result of this normalisation process is to provide normalised x- 
component traces that are approximately consistent across the receiver spread, 
and that may be compared directly to one another. This process effectively 
enables the x-component of the particle iiibtion recorded "at areceiver to be 
calibrated against the vertical component of the particle motion recorded at that 
receiver (which should contain predonainantly p-wave events and so be 
relatively unaffected by the near-surface or coupling variations). 
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If the y-component of the particle motion at the receiver is acquired as well as, 
or in addition to, the x-conaponent of the particle motion, it would be possible to 
calibrate the y-component of the particle motion against the vertical component 
of the particle motion. 

If full 3 C or 4C seismic data are available, it is possible to compute two cross- 
correlograms - a first cross-correlogram of the x-component against the z- 
- component of the particle motion, and a second cross-correlogram of the y- 
component against the z-component of the particle motion. Each of these 
correlograms may then be normalised so that the peak in the correlogram at zero 
time delay has the same amplitude in each trace. The two normaUsed cross- 
correlograms may then be directly compared against each other - because this 
process calibrates the x-component of the particle motion against the z- 
component of the particle motion and also calibrates tiie y-component of the 
particle motion against the z-component of the particle motion, and so is 
equivalent to calibratmg.the x-component of the.particle motion against the y- 
component of the particle motion. 

Figure 7 includes only information about the amplitude variations of the traces 
in Figure 6. However, it can be seen that the waveforai of the "statics peak" in 
the cross-correlogram traces of Figures 6(b) varies from one receiver position to 
another. The variations in the waveform of the static peaks may also provide 
information about the near surface and/or variations in coupling strength 
between the receiver and the earth from one receiver position to another. 

Conceptually, it would be possible to construct a figure showing the variation in 
some measure of flie waveform of the "statics peak*' - such as, for example, the 
the fiill width at half maximum height of the wavefi)rm - in the normalised 
receiver functions of Figures 6(b), Such a figure would show how the waveform 
of this peak varies from one receiver position to another, and these variations hi 
the waveform of the "statics peak" could also provide information about the 
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geological effect of the near-surface and about variations in coupling strength 
between the receiver and the earth. 

Alternatively, it would be possible to derive a jSlter that transforms the 
wavefonn of the "statics peak** in one of the receiver functions of Figure 6(b) to 
an undistorted waveform, such as the wavefonn of the peak at zero time lag in 
the trace. The filter may then be applied to the original seismic data, for 
example to the original trace of the x-component of the particle motion. This 
process may be repeated for each receiver - the trace of the x-component for 
each receiver is transformed by the filter derived from the wavefonn of the 
"statics peak" in the receiver fimctions for that receiver. This process should 
provide normalised x-component traces that are approximately consistent across 
the receiver spread, and that may be compared directly to one another. 

In the embodiments described above the receiver fimctions have been calculated by 
cross-correlation. The invention is not hmited to this, however, and any suitable 
correlation technique or algorithm may be used to determine the receiver fimctions. For 
example the receiver fimctions could be obtained by deconvolution between component, 
and the amplitude and waveform information in flie statics peak can then still be used in 
the same way as described above. Indeed, a receiver fimction obtained by 
deconvolution may be preferable since a deconvolution effectively removes the specific 
source signature of the p-wave (as it exists upon incidence at the near-surface/basement 
interface) such that the only effects left in the receiver fimction are those due to near- 
surface structure and/or vector-infidelity. 

A fiirther example of a suitable technique for obtaining the receiver fimctions is a 
bicoEerence tirhe"delay esfimatia as described by X. Ikelfe in "Geophysics", 

Vol62,pl947 (1997). 



In the description of the invention given above, components of the particle 
velocity acquired at the receiver have been used as the components of particle 
motion in calculating the receiver fimctions. In principle, however, components 



MAC Folio No P51798GB ' 

of the particle displacement at the receiver, or coinponents of the particle 
acceleration at the receiver, or higher derivatives of the particle velodty, could 
alternatively be used. 

Indeed, the invention is not limited to use of components of the particle motion, 
and any suitable parameters of the acquired seismic data may be used to obtain 
the receiver functions. For example, the acquired seismic data may be 
deconq>osed into its up-going and down-going components, or into its p- and s- 
components, and these may be used to obtain the receiver functions. 

The invention may also be ^plied to p- and s-wave "reverberation events". In 
such events partial mode conversion, for example upon transmission through or 
reflection at the interfece between the near surface 5 and the basement 6. In a 
reveibraation event, however, the p- and s-waves do not pass directly fix>m the 
point at which the mode-conversion occurs to the receiver. In a reverberation 
.eventj-however, the p.- and^waves imdergo one or.more.reflection at a 
boundary of the near-surface and so make multiple passes through the near- 
surface 5 before being incident on the receiver. 

Seismic data may also contains events that arise from partial mode conversion 
that occurs not at a boundary of the near-surface but within the near-surface. 
This may arise owing to, for example, layering effects in the near-surface. 

In the embodiments described above one of the p- and s- wave events has been 
generated by partial mode conversion. However, as noted above, the invention is not 
limited to this and may be ^lied to any pair of corresponding p- and s-events for 
which differences in amplitude and/or waveform of the two events arise primarily from 
the near-surface 5. For example, in a seismic survey that uses a source that emits both 
p- and s-waves, the acquired data will contain an event arising from reflection of p- 
waves at a particular point on the lower boundary of the near-surface 5 and will also 
contain a corresponding event arising from reflection of s-waves at that point on the 
lower boundary of the near-surface 5. The invention may be appUed to these events. (It 
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should be noted that, since a water column will not support s-wave propagation, a 
source emitting both s- and p-modes would need to be disposed on the seabed, on the 
earth's surface or in a borehole. 

Figure 8 is a block flow diagram illustrating one embodiment of the present 
invention. Initially, at step 1, multi-component seismic data is acquired at a 
survey location. The seismic data should include the vertical component of the 
particle motion and, at least one horizontal component of the particle motion. 
The data is acquired at a plurality of receivers disposed at diflFerent locations in 
the survey area. The multi-component data may be fiill 3C or 4-C seismic data. 

Rather than applying the invention to newly-acquired seismic data, the invention 
may also be applied to pre-existing survey data. Step 1 may therefore be 
replaced by the step, step 2, of retrieving suitable multi-component seismic data 
from storage. 

At step 3 seismic data relating to the i*^ receiver is selected. At step 4 the 
receiver ftmction for the i* receiver is calculated by from two parameters of the 
seismic data acquired at the i* receiver. In this embodiment the receiver 
ftmction is found by cross-correlating flie vertical component of the particle 
" motion acquired at the i* receiver is cross-correlated with a horizontal 
component of the particle motion acquired at the i* receiver. In tihis 
embodiment it is assumed that the horizontal component is the x-component 
(radial component) of the particle motion at the receiver and that the particle 
motion recorded at the receiver was the particle velocity, but the invention is not 
limited to this. 



The results of step 4 give a receiver ftmction of the type shown in Figure 5 or 
6(a), for the i^ receiver. 
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At Step 5 the receiver function for the i* recdver is normalised, relative to the 
amplitude of the peak in the receiver function at zero time lag. The result of step 
5 is the normalised receiver function as shown in Figure 6(b) for the i* receiver. 

At stq) 6 the "statics peak" in the normalised receiver function is determined. 
This may be done nsing an automatic picking method, for example, or it may be 
done by outputting the results of step 5 for display by an operator who manually 
picks the "statics peak". 

At step 7 the ampUtude of the "statics peak" in the normaHsed receiver function 
is determined. At step 8, the ampUtude of the "statics peak" is used to normaHse 
the origiaal trace of the horizontal component of the particle motion recorded at 
the i* receiver. 

At step 9 it is determined whether the data ftom all receivers in the receiver 
array has been processed. If the result of the determination is "no", a suitable 
counter is increased by one at step 10, and steps 3 to 9 are then repeated for the 
(i+l)**^ receiver. 

When a "yes" determination is obtained at step 9, this indicates that data from all 
receivers in the array has been processed. The results of the processing may 
then be output at step 1 1 . The data may be output in any suitable form - for 
example it may be printed onto p^er, displayed on a visual display unit, or it 
may be stored in any suitable form for subsequent retrieval. Alternatively, the 
data may be subjected to additional processing to determine information relating 
to the geological structure of the survey location. 

The foim in which data is output at step 11 can be chosen to suit the particular 
needs of the user. The data that is ou^ut may include some or all o'f the receiver 
function produced at step 4, the normalised receiver function produced at step 5, 
the amplitude of the "static peak" in the normalised receiver function as 
deteraiined at step 7, or the normalised trace of the horizontal component of the 
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particle motion as detennined at step 8. These data may be output for each 
receiver, or only for selected receivers. 

The method of the invention is not limited to performing the steps in the order 
shown in Figure 8. For example, step 6 of picking the "statics peak" can be 
carried out before step 5 of normalising the receiver function. Furthermore, in 
the method shown in Figure 8 steps 3 to 8 are performed in sequence for one 
receiver, and are then repeated for another receiver and so on. It would 
alternatively be possible to perform steps 3 and 4 for one receiver, and then 
repeat steps 3 and 4 for other receivers before moving onto step 5, Similarly, 
step 5 of normalising the receiver function could be performed for each receiver 
before moving onto step 6 of picking the "statics peak". Likewise, step 6 of 
picking the "statics peak*' could be performed for each receiver before moving 
on to step 7 of detennining the amplitude of the "statics peak'*, and step 7 of 
determining the amplitude of the "statics peak*' could be performed for each 
receiver before moving onto step 8 of normalising the horizontal component of 
the particle motion. 

The method of Figure 8 is suitable if it is possible to identify the "statics pealsf' or 
another event relating to the traveltime difference between p- and s-waves through the 
near-surface in flie individual, unstacked receiver functions. In some cases, however, it 
will not be possible to identify the "statics peak" in the individual receiver functions 
and, in such cases, it is necessary to include a step of stacking receiver functions to 
improve the signal-to-noise ratio. 

If a stacking step is necessary to allow identification of the "statics peak", step 4 of 
Figure 8 would be replaced by the stqp of calculating all the receiver functions for the i* 
receiver gather. The i* receiver gather contains all traces recorded at the i* receiver and 
for each of tihiese traces it is possible to calculate a receiver function. This step is 
followed by the step (not shown in Figure 8) of stacking all receive functions that have 
been calculated for the i**^ receiver gather. Steps 5 to 7 of Figure 8 are then carried out 
on the stacked receiver function for the i* receiver gather. 
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The stacking process is then repeated for the receiver functions determined for other 
receiver gathers. 

Figure 9 is a schematic block diagram of a programmable apparatus 15 
according to the present mvention. The system comprises a programmable data 
processor 8 with a program memory 9, for instance in the form of a read only 
memory ROM, storing a program for controlling the data processor 8 to _ 
perform, for example, a method as illustrated in Figure 8. The system further 
comprises non-volatile read/write memory 10 for storing, for example, any data 
which must be retained in the absence of power supply. A * Vorking" or 
"scratchpad" memory for the data processor is provided by a random access 
memory (RAM) 11. An input interface 12 is provided, for instance for receiving 
commands and data. An output interface 13 is provided, for instance for 
displaying information relating to the progress and result of the method. User- 
dejBned data may be supplied via the mput interface 12 or may optionally be 
provided by a machine-readable store 14. 

The program for operating the system and for performing the method described 
hereinbefore is stored in the program memory 9, which may be embodied as a 
semi-conductor memory, for instance of the well-known ROM type. However, 
the program may be stored in any other suitable storage medium, such as 
magnetic data carrier 9a (such as a "floppy disc") or CD-ROM 9b. 
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CLAIMS: 

1. A method of processing seismic data including corresponding first and second 
modes of seismic energy, the method comprising the steps of: 

a)/ processing a trace indicative of a jSrst parameter of tiie seismic data acquired at a 
jSrst receiver including at least events corresponding to the first mode with a trace 
indicative of a second parameter of the seismic data acquired at the first receiver 
including at least events corresponding to the second mode; 
h) identifying an event in the processed data corresponding to partial mode 
conversion; and 

c) obtaining information firom the amplitude and/or the waveform of the event in 
the processed data corresponding to partial mode conversion. 

2. A mefliod as claimed in claim 1 wherein the second mode has been generated by 
partial mode conversion of the first mode at a boundary face of a layer of the earth. 

3. A method as claimed in claim 1 or 2 wherein step (c) comprises normalising the 
amplitude in the processed data of the event corresponding to partial mode conversion, 

4. A method as claimed in claim 3 wherein step (c) comprises normalising the 
amplitude in the processed data of the event corresponding to partial mode conversion 
relative to the ampUtude of the processed data at zero time delay. 

5. A method as claimed in claim 3 or 4 and comprising the fiarther step of 
normalising flie trace indicative of tiiie second parameter relative to the normalised 

"amplitude of the of the event corresponding to partiail mode conversion: " 

6. A method as claimed in claim 1, 2, 3, 4 or 5 and comprising the further step of 
rq)eating steps (a), (b) and (c) for seismic data acquired at a second receiver, the second 
receiver not being co-located with the first receiver. 
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7. A mefliod as claimed in any preceding claim wherein the first parameter of the 
seismic data is a vertical component of particle motion acquired at the respective 
receiver. 

8. A mettiod as claimed in claim 7 wherein the second parameter of the seismic 
data is a horizontal component of particle motion acquired at the respective receiver 

9. A method as claimed in claim 7 wherein the second parameter of the seismic 
data is the component of particle motion acquired at the respective receiver in the 
direction between a source of seismic energy and the respective receiver. 

10. A method as claimed in any of claims 1 to 6 wherein at least one of the first and 
second parameters of the seismic data is obtained by decomposing the seismic data 
acquired at the respective receiver. 

1 L A method as claimed in any preceding claim wherein step (a) comprises cross- 
correlating the trace indicative of the first parameter of the seismic data acquired at the 
first receiver with the trace indicative of the second parameter of the seismic data 
acquired at the first receiver, 

12. A method as claimed in any of claims 1 to 10 wherein step (a) comprises 
deconvolving the trace indicative of the first parameter of the seismic data acquired at 
the first receiver from the trace indicative of the second parameter of the seismic data 
acquired at the first receiver. 

13. A method of seismic surveying comprising flie stqps of: directing seismic energy 
propagating in a first mode towards a boundary face of a layer of the seabed such tha.t 
partial mode conversion of the seismic energy into a second mode occurs at the 
boundary face; acquiring seismic data including the first and second modes of seismic 
energy at one or more receivers; and processing the seismic data according to a method 
defined in any of claims 1 to 12. 
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14. An apparatus for processing seismic data including corresponding jSrst and 
second modes of seismic energy, the apparatus comprising: 

a) means for processing a trace indicative of a first parameter of the seismic data 
acquired at a first receiver including at least events coiresponding to the Sist mode with 
a trace mdicative of a second parameter of the seisndc data acqxiired at the first receiver 
including at least events corresponding to the second mode; 

b) means for identifying an event in the processed data corresponding to partial 
mode conversion; and ' 

c) means for obtaining information firom the ampUtude and/or the waveform of the 
event in the processed data corresponding to partial mode conversion, 

15. An apparatus as claimed in claim 14 and comprising normalising means for 
noimalising the amplitude in the processed data of the event corresponding to partial 
mode conversion. 

16. An apparatus as claimed in claim 15 wherein the normalising means are adapted 
to normalise the amplitude in the cross-correlated data of the event corresponding to 
partial mode conversion relative to the amplitude of the cross-correlated data at zero 
time delay. 

17. An apparatus as claimed in any of claims 14 to 16 and comprising means for 
normalising the trace indicative of the second parameter relative to the normalised 
amphtude of the of the event corresponding to partial mode conversion. 

18. An apparatus as claimed in any of claims 14 to 17 and comprising a 
programmable data processor. 



19. A storage medium containing a program for a data processor of an apparatus as 
claimed ia claim 18. 
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ABSTRACT 
A method of processing seismic data 



A method of processing seismic data that includes including first and second modes of 
seismic energy where the second mode has been generated by partial mode conversion 
of the first mode at aboundary face of a layer of the seabed includes the step of cross- 
correlating a trace acquired at a receiver and including events corresponding to the first 
mode with a trace acquired at "the same receiver and including events corresponding to 
the second mode. An event in the cross-correlated data that corresponds to partial mode 
conversion is identified, and the ampUtude of this event is normahsed, for example 
relative to the amplitude of the peak in the cross-correlogram at zero time delay. 

Information about the effects of the static shift produced by the layer and/or about 
vector infideUty can be arrived firom normalised cross-correlograms for receivers in a 
receiver array. 




TH/S PAGE BUNK 



(usno) 



THIS PAGE BUNK (USPTO) 




Up«gomg F 



, Schematic of noiseless traces recorded for the situation in Figure 2 where a 
single up-coming P-arrival transmits and converts to an S-anival at the ^ur(<tce layer 
interface. The time difference between principal P- and S-arrivals is dt The XKJomponent 
has been amplified. 
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Cross-correlogram between vertical and horizontal components in Figure?. The 
peak at time 0 is due to leakage between components as shown in Figure3v whereas the 
peak at time dt is due to the P-S conversion. 

Leakage Conversion 
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Robertsson and Curtis (2000). 
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Receiver nr. 



Profile of full CRO P*VX crosscorrelatfon stack (normalized) 




Receiver 



Figures 6(a) and 6(b): Same data as shown in Figure 5- Figure 6(a) 
normalizing to the peak at lag-time 0. 



(Upper): raw data. Figure 6(b) (Lower): data after 
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Figure 7: A perturbation correction estimated from the plots in Figure 6(b). 
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